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Abstract. We analyze a spatially homogeneous SU(2) Yang-Mills-Higgs 
system both in classical and quantum mechanics. By using the Toda crite- 
rion of the Gaussian curvature we find a classical chaos-order transition as 
a function of the Higgs vacuum, the Yang-Mills coupling constant and the 
energy of the system. Then, we study the nearest-neighbour spacing dis- 
tribution of the energy levels, which shows a Wigner-Poisson transition by 
increasing the value of the Higgs field in the vacuum. This transition is a 
clear quantum signature of the classical chaos-order transition of the system. 



1 Introduction 



It is well known that the spatially uniform limit of scalar electrodjTiamics 
and Yang-Mills theory exhibits classical chaotic motion [1-8]. Usually the 
order-chaos transition in these systems has been studied numerically with 
Lyapunov exponents and sections of Poincare. Less attention has been paid 
to analytical criteria. 
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In this work we study analytically the suppression of classical chaos in 
the spatially homogenous SU(2) Yang-Mills-Higgs (YMH) system induced 
by the Higgs field. Then we analyze the energy fiuctuation properties of the 
system, which give a clear quantum signature of the classical chaos-order 
transition of the system [9-13]. 

The SU(2) YMH system describes the interaction between a scalar Higgs 
field (f) and three non-Abelian Yang-Mills fields A^, a = 1,2,3. The La- 
grangian density of the YMH system [14] is given by 

L = l{D,<p)^D^<p) - v{<p) - , (1) 

where 

(2) 

F;, = d.At - d.A-^ + ge-'^AlAl , (3) 

with = (T^/2, b = 1,2,3, generators of the SU(2) algebra, and where the 
potential of the scalar field (the Higgs field) is 

V{cj>) = i^'lcPf + X\cl>\\ (4) 

We work in the (2+l)-dimensional Minkowski space (/i = 0, 1, 2) and choose 
spatially homogeneous Yang-Mills and the Higgs fields 

diAl^di<f>^0, i = l,2 (5) 

i.e. we consider the system in the region in which space fiuctuations of fields 
are negligible compared to their time fiuctuations. 

In the gauge Aq — and using the real triplet representation for the Higgs 
field we obtain 

L = ^ + ^{A, + A,) - g'^-AUl - \{A, . A,f+ 

+ + - (A • - {A, . $f\ - V($) , (6) 

where = (f)^ , f^) , Ai = {A\,Al,Al) and A2 = {AlAlAl). 

When //^ > the potential V has a minimum at |0[ = 0, but for /x^ < 
the minimum is at 
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which is the non zero Higgs vacuum. This vacuum is degenerate and after 
spontaneous symmetry breaking the physical vacuum can be chosen 0o = 
(0,0,^;). If A\ — qi, Al — 52 and the other components of the Yang-Mills 
fields are zero, in the Higgs vacuum the Hamiltonian of the system reads 

H-lipi+ pI) + g'v^l + ql) + lo'qlql , (7) 

where pi — qi and p2 — q2- Here — 2g^v^ is the mass term of the 
Yang-Mills fields. This YMH Hamiltonian is a toy model for classical non- 
linear dynamics, with the attractive feature that the model emerges from 
particle physics. In the next sections we analyze first the classical chaos- 
order transition of the YMH system and then its connection to the quantal 
fluctuations of the energy levels. 



2 Classical transition from chaos to order 

In this paper we study the chaotic behaviour of this YMH system by using 
the Toda criterion of the Gaussian curvature of the potential energy [16]. 
The Toda criterion is based on a local estimation of the rate of separation of 
neighboring trajectories in the classical phase space of the model. To obtain 
the time evolution of a dynamical system with the Hamiltonian 

H-§ + § + V{q„q2), (8) 



the following Hamilton equations have to be solved: 

dcL _dH dp _ _dH 
dt dp ^ dt 9q ' 



(9) 



where q = (51,52) and p = (pi,P2)- The linearized equation of motion for 
the deviations are 

f = /.p, f = (10, 

where / is the 2x2 identity matrix, and 

^''^^^ = %^'''='^^*^ ' ^^^^ 
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with q(t) sohition of the Hamihon equations. The stabihty of the dynamical 
system is then determined by the eigenvalues of the 4x4 stability matrix 



r(qW) 




(12) 



If at least one of the eigenvalues of the stability matrix T is real, then the 
separation of the trajectories grows exponentially and the motion is unstable. 
Imaginary eigenvalues correspond to stable motion. 

To diagonalize the matrix F, we must first solve the equation of motion. 
The problem can be significantly simplified by assuming that the time depen- 
dence can be eliminated by replacement of the time-dependent point q(t) of 
configuration space by a time-independent coordinate q, i.e. F(q(t)) = F(q). 
The eigenvalues then are 



A 



where 



B 



+ 



C 



dql 



dql 



(13) 



(14) 



:i5) 



dql dql ^dqiq2' 

Now, if > then with C > the eigenvalues arc purely imaginary and the 
motion is stable, while with C < the pair of eigenvalues becomes real and 
this leads to exponential separation of neighboring trajectories, i.e. chaotic 
motion. The parameter c has the same sign as the Gaussian curvature Kq 
of the potential-energy surface: 



Kg{(1u(12) 



dq'i dq'^ 



dq\qi ' 



For our YMH system the potential energy is given by 



^(?i, 52) = ^'^^'(gi + (&) + 2/^1^2 



and the F matrix reads 



r(gi,g2) 












-25'^gig2 

-2^2y2 _ ^2^2 



\ 
1 



/ 



(16) 



(17) 



(18) 
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At low energy, the motion near the minimum of the potential, where the 
Gaussian curvature is positive, is periodic or quasi-periodic and is sepa- 
rated from the instability region by a line of zero curvature; if the energy 
is increased, the system will be, for some initial conditions, in a region of 
negative curvature, where the motion is chaotic. According to this scenario, 
the energy E^. of chaos-order transition is equal to the minimum value of the 
line of zero gaussian curvature K{qi, ^2) on the potential-energy surface. For 
our potential the gaussian curvature vanishes at the points that satisfy the 
equation 

{2g\^ + g\l){2g\^ + g^^) - Ag^lql = . (19) 

It is easy to show that the minimal energy on the zero-curvature line is given 
by: 

Ec = VrniniKo = 0, gi) = Qg'v^ , (20) 

and by inverting this equation we obtain Vc = {E /Gg'^Y^^. We conclude that 
there is a order-chaos transition by increasing the energy E of the system 
and a chaos-order transition by increasing the value v of the Higgs field in 
the vacuum (see also [2]). Thus, there is only one transition regulated by the 
unique parameter E/{g'^iA). 

It is important to point out that in general the curvature criterion guar- 
antees only a local instability [16] and should therefore be combined with 
the Poincare sections [17-18]. In the paper [19] we used a fourth-order 
Runge-Kutta method to compute numerically the classical trajectories and 
the Poincare sections. The results confirm the analytical predictions of the 
curvature criterion: with E = 10 and = 1 we get the critical value of the 
onset of chaos Vc = {E/Qg'^Y^'^ ~ 1.14. 

3 Quantum chaos and the Wigner distribu- 
tion 

The energy fluctuation properties of systems with underlying classical chaotic 
behaviour and time-reversal symmetry agree with the predictions of the 
Gaussian Orthogonal Ensemble (GOE) of random matrix theory, whereas 
quantum analogs of classically integrable systems display the characteristics 
of the Poisson statistics [9-12]. Some results in this direction for field theo- 
ries have been obtained by Halasz and Verbaarschot: they studied the QCD 
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lattice spectra for staggered fermions and its connection to random matrix 
theory [13]. 

In quantum mechanics the generahzed coordinates of the YMH system 
satisfy the usual commutation rules = iSki, with k,l — 1, 2. Introduc- 

ing the creation and destruction operators 





H ^ \ nQk - i\ 7r-Pk : (21) 



the quantum YMH Hamiltonian can be written [15] 

^ = ^0 + Ig'v , (22) 

where 

Hq = u;{afai + aja2 + 1) , (23) 

V^^{ai + atf{a2 + atf, (24) 

with u;^ = 2g'^v'^ and [at, af] = 6ki, k,l = 1,2. 

The most used quantity to study the local fluctuations of the energy 
levels is the spectral statistics P{s). P{s) is the distribution of nearest- 
neighbour spacings Sj = {Ei+i — Ei) of the unfolded levels Ei. It is obtained 
by accumulating the number of spacings that lie within the bin (s, s + As) 
and then normalizing P{s) to unit [9-12]. 

For quantum systems whose classical analogs are integrable, P{s) is ex- 
pected to follow the Poisson hmit, i.e. P{s) = exp (— s). On the other hand, 
quantal analogs of chaotic systems exhibit the spectral properties of GOE 
with P{s) = (7r/2)s exp (— |s^), which is the so-called Wigner distribution 
[9-12]. The distribution P{s) is the best spectral statistics to analyze shorter 
series of energy levels and the intermediate regions between order and chaos. 

Seligman, Verbaarschot and Zirnbauer [20] analyzed a class of two-dimensional 
anharmonic oscillators with polynomial perturbation by using the Brody dis- 
tribution [21] 

P{s,u;) = a{uj + l)s'^exp{-as'^+^) , (25) 



with 



a-{n^]r'- (26) 
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This distribution interpolates between the Poisson distribution {uj = 0) of 
integrable systems and the Wigner distribution {u = 1) of chaotic ones, and 
thus the parameter cu can be used as a simple quantitative measure of the 
degree of chaoticity. 



Higgs vacuum v 


Brody parameter uo 


1.00 


0.99 


1.05 


0.47 


1.10 


0.34 


1.15 


0.12 


1.20 


0.01 



The Table shows the calculated Brody parameter co of the P{s) distri- 
bution for different values of the Higgs vacuum v. There is Wigner-Poisson 
transition by increasing the value v of the Higgs field in the vacuum. Thus, 
by using the P(s) distribution and the Brody function, it is possible to give a 
quantitative measure of the degree of quantal chaoticity of the system. Our 
numerical calculations show clearly the quantum chaos-order transition and 
its connection to the classical one. 

We compute the energy levels with a numerical diagonalization of the 
truncated matrix of the quantum YMH Hamiltonian in the basis of the har- 
monic oscillators [22]. If \n1n2 > is the basis of the occupation numbers of 
the two harmonic oscillators, the matrix elements are 

< r^^r^2|^o|r^lr^2 >= w(ni + + 1)5„;„^5„^„, , (27) 

and 

< n[n^\V\nin2 >= ^[^ni{ni - 1)5„;„^_2 + + + 2)5„;„^+2+ 

+ (2ni+l)5„;„Jx[^n2(n2 - ^(^2 + l)(n2 + 2)5„;„^+2+(2^2+l)5„;„J ■ 

(28) 

The symmetry of the potential enables us to split the Hamiltonian matrix 
into 4 sub-matrices reducing the computer storage required. These sub- 
matrices are related to the parity of the two occupation numbers rii and 
n2'- even-even, odd-odd, even-odd, odd-even. The numerical energy levels 
depend on the dimension of the truncated matrix: we compute the numerical 
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levels in double precision increasing the matrix dimension until the first 100 
levels converge within 8 digits (matrix dimension 1156 x 1156) [22-23]. 

We use the first 100 energy levels of the 4 sub-matrices to calculate the 
P{s) distribution. In order to remove the secular variation of the level density 
as a function of the energy E, for each value of the coupling constant the 
corresponding spectrum is mapped, by a numerical procedure described in 
[24], into one which has a constant level density. 

4 Conclusions 

The chaotic behaviour of an homogenous YMH system has been studied both 
in classical and quantum mechanics. The Gaussian curvature criterion shows 
that the chaotic behaviour is regulated by the unique parameter E/{g'^v^). 
The YMH system has a order-chaos transition by increasing the energy E 
and a chaos-order transition by increasing the value v of the Higgs field in 
the vacuum. 

The -P(s) distribution of the energy levels confirms with great accuracy 
the classical chaos-order transition of the YMH system. In particular, the 
Brody parameter uo shows a Wigner-Poisson transition for the P{s) distri- 
bution in correspondence to the classical chaos-order transition. 
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